Quantum Mechanical Limits to Inertial Mass Sensing by Nanomechanical Systems 
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We determine the quantum mechanical limits to inertial mass-sensing based on nanomechanical 
systems. We first consider a harmonically oscillating cantilever whose vibration frequency is changed 
by mass accretion at its surface. We show that its zero-point fluctuations limit the mass sensitivity, 
for attainable parameters, to about an electron mass. In contrast to the case of a classical cantilever, 
we find the mass sensitivity of the quantum mechanical cantilever is independent of its resonant 
frequency in a certain parameter regime at low temperatures. We then consider an optomechanical 
setup in which the cantilever is reflective and forms one end of a laser-driven Fabry-Perot cavity. For 
a resonator finesse of 5 the mass sensitivity at T = is limited by cavity noise to about a quarter 
of a Dalton, but this setup has a more favorable temperature dependency at finite temperature, 
compared to the free cantilever. 

PACS numbers: 07.10.Cm, 42.50 Lc, 42.50.Pq 



Mechanical oscillators are quickly approaching the 
quantum regime. Advances in nanofabrication, laser- 
cooling, electromechanics etc. have allowed experimen- 
talists to cool harmonically oscillating cantilevers to an 
excitation number of n ~ 25 the eventual aim be- 
ing to reach the ground state n = 0, see Ref. Q and 
references therein. The availability of macrosocopic me- 
chanical systems operating in the quantum regime will 
allow one to address conceptual questions regarding the 
nature of the quantum-classical interface [3(, to coher- 
ently control atomic [H and molecular [|| systems; and 
to construct displacement [f| , mass Q , and force Q sen- 
sors of unprecedented sensitivity. 

A prime example of an ultra-precise displacement sen- 
sor is the Laser Interferometer Gravitation- Wave Obser- 
vatory (LIGO), which is expected to enter the quantum 
regime in the coming years A lar ge body of work, 
both experimental [9( and theoretical [10| has been de- 
voted to investigating the ultimate, quantum mechanical 
limits to the displacement sensitivity of devices such as 
LIGO. 

Mass sensing is another important application of 
nanomechanical resonators useful in physical, chemical 
and biological contexts 0]. Typically the mass to be 
"weighed" SM is deposited on a harmonically oscillating 
cantilever of mass M SM , whose resonant frequency 
wo changes as a result. By successfully detecting that 
small change Su>o, sensitivities have been reached where 
a single molecule can be weighed 

The question that we address in this Letter is to de- 
termine the ultimate resolution of nanocantilever-based 
mass sensing. As is the case for gravitational wave de- 
tectors this limit is set by the laws of quantum mechan- 
ics. However, to the best of our knowledge all previous 
descriptions of cantilever-based mass-sensing have been 
classical 
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or have accounted for noise 
effects only at high temperature |l5J]. These treatments 
predict a sensitivity to infinitesimally small masses as the 



temperature is lowered (T — > 0), but of course, this result 
has to be corrected to account for the effects of quantum 
noise. Motivated by the rapid approach of mechanical 
oscillators to the quantum regime, we present a simple 
quantum mechanical analysis of mass sensing illustrating 
the limits set by thermal noise at low temperatures. Our 
results confirm that in the T — > limit, mass sensitivity 
is limited by purely quantum mechanical effects. 

We begin by considering the simple case of a mechani- 
cal cantilever of mass M and vibration frequency ujq that 
is driven near resonance by an electrical circuit, typically 
a phase-locked loop that can detect its frequency shift 
8loq when a mass SM is deposited on it (see Ref. [7| for a 
diagram of the circuit). For mass increments SM <C M 
we have Ml 



oM = - — duj , 
duo 



(1) 



where the mass responsivity dM/duo = —2M/u) 

Various noise sources also couple to the cantilever, re- 
sulting in random frequency shifts. In order to obtain a 
signal to noise ratio (SNR) of at least unity the smallest 
detectable frequency shift Sojq iS should equal the mean 
frequency shift due to the various noise sources, 



Su)Q. s = 



Su(u})T(w)du) 
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where S u (uj) is the spectral density of noise-induced fre- 
quency fluctuations and T(uj) is the normalized transfer 
function of the circuit realizing the measurement. Follow- 
ing , that expression can be simplified if the measure- 
ment circuit responds only to frequencies in a bandwidth 
A/ around loq. Approximating then T(u>) by an ideal 
band pass filter with T(oS) = 1 between loq — ttA/ and 
ljq + 7rA/ and elsewhere, Eq. ([2|) yields the smallest 
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detectable frequency shift as 



/ S u {uj)d,uj 

Juo—irAf 
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The spectral density Su(uj) depends on the specific sys- 
tem noise source being considered. In order to isolate the 
fundamental noise limitations to mass sensing we assume 
in the following that both the drive and the detection 
arc noiseless, despite the fact that in reality they intro- 
duce backaction noise. In practice, transduction needs to 
be optimized to reduce backaction noise to a level com- 
parable to the intrinsic noise that we consider here (l4j |. 
Considering the parameters given in the caption of Fig. [TJ 
we also neglect other technical sources of noise such as 
temperature fluctuations, which yield limitations to mass 
sensitivity smaller than one electron mass for tempera- 
tures below IK 0; elastic and inelastic collisions with 
background gas particles, whose influence on mass sensi- 
tivity also drops below one electron mass for a vacuum of 
10~ 6 Torr an d defect motion within the cantilever 

itself. As pointed out in Ref. [3] , under these conditions 
the operational limits of mass sensing are set by thermo- 
mechanical displacement fluctuations in the cantilever. 

With these caveats in mind we treat the cantilever as 
a quantum harmonic oscillator of mechanical quality fac- 
tor Q and total thermal noise energy E no i Se = 5 K (<7noisc)j 
where k = Mu>q is the stiffness of the cantilever and 
(^noise) ^ s ^s center-of-mass mean square displacement. 
The cantilever is driven on resonance by a noiseless cir- 
cuit with driving energy Ed- 

From the equation of motion of the cantilever, the spec- 
tral density of noise-driven random displacements S q (oj) 
is given by [13] 

2u)Q -Encase 
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(4) 



This spectral density can be linked to the spectral 
density of frequency fluctuations S^(ui) through the 
relations S q (uj) = (q%) S '^{lo) / '(2it) 2 and ~ 
uj 2 S c i,(uj)/(2tt) 2 Q 2 , where S^,(lu) is the spectral density 
of phase fluctuations and (g 2 ,) is the mean squared am- 
plitude of the cantilever resulting from the (noiseless) 
resonant drive jl8j |. Combining these expressions with 
Eqs. ([3]) and ([4J the smallest mass detectable is given by 



where n — ( e huj o/k B T _ j s ^ nc mcan number of ther- 
mal phonons at frequency tuo and temperature T. E no \ se 
can be easily modified to account for cantilever damping, 
see e.g. Ref. [l9[. However these corrections are negligi- 
ble for the high quality factors being considered here for 
mass sensing. At T = Eq. ([6]) gives E no [ se = fuvo/2, 
the ground state energy of the cantilever, while at high 
temperatures T 3> hw^/kB, E noisc = kgT which yields 
the classical results published earlier. 

The integral from Eq.© has been evaluated analyti- 
cally in the limit Q ^ 1 and plotted in Fig. [TJ An illumi- 
nating form for this integral is obtained with the further 
approximation 2irAf <C i^o/Q (this can be achieved us- 
ing the parameters of Fig. Q] but Q = 10 5 ) simplifying 
Eq.© to 
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where Ed = ^A/u^g 2 ,) is the driving energy. 

To evaluate 5M S explicitly requires an expression for 
the noise energy -Excise ■ For a quantum harmonic oscil- 
lator in equilibrium with a thermal reservoir at temper- 
ature T we have 

E noisc = Tiw {n + i), (6) 



S1VL = 2V2irM 



huJo(2n + 1) 
2E d 
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As expected, quantum fluctuations prevent the smallest 
detectable mass 5M S from going to zero when T — ► 0. 
Rather, it is limited by the quantum zero-point fluctua- 
tions of the cantilever to 



5M S (T = 0) = 2y/^M (J^J 
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Quj 
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which yields 8M S (T ~ 0) ~ 10 electron masses for our 
parameters. Note that we have retained ojq in Eq. ([5]) in- 
side the first parentheses to provide a scale for the driving 
energy, and inside the second parentheses to facilitate 
comparison with the corresponding classical expression 
0]. However these frequency dependencies cancel each 
other in the quantum regime, indicating that in contrast 
to the classical regime mass sensitivity is now indepen- 
dent of the mechanical frequency of the cantilever within 
the range of validity of our description. 

Examples of quantum and classical mass sensitivities 
are plotted as a function of temperature in Fig. [T] For 
the parameters of that figure, we have SAI S (T = 0) ~ 
5 x 10~ 31 kg, or about an electron mass. This limit is 
reached when n <C 1, which occurs for ojq ~ 1GHz 
at T < lOmK. Such temperatures have recently been 
reached experimentally for mechanical oscillators using 
laser cooling techniques [iFj| . These same techniques also 
hold the promise of mechanical squeezing of the can- 
tilever motion that would result in E nolsc < hiuo/2 (21. 22| 
with an associated improvement in 6M S . 

We now turn to an analysis of the mass sensitivity 
afforded by a typical optomechanical cooling configura- 
tion where the vibrating cantilever forms one end of a 
Fabry-Perot cavity driven by laser radiation. The can- 
tilever itself is driven as before by an external noiseless 
circuit with energy Ed- Let p and q be respectively the 
momentum and displacement operators of the cantilever 
with [q,p] = ih. In terms of the dimcnsionlcss variables 
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FIG. 1: (Color online) Quantum (solid curve) and classical 
(dashed curve) mass sensitivity 8M a in electron mass units 
(1EMU ~9x l(T 31 kg) of the free cantilever mass sensor for 
M ~ O.lfg, Q = 10 7 , u ~ 1GHz, A/ ~ 1kHz and E d ~ 
1.6 x KT 15 J. 



p = pj \/HMu! and q = q^J AIuj /h, the Hamiltonian that 
models that cavity mass sensor is 



via) a + (p 2 + <f ) - hga)aq - E d (t)q. (9) 



where a and a) are bosonic annihilation and creation 
operators for the cavity mode of frequency u>i, luq is 
the natural vibration frequency of the cantilever and 
g = (uji/L ) (h/ Mujr,) 1 / 2 is the optomechanical coupling 
parameter [23j. 

To find the displacement fluctuations of the micro- 
mirror due to noise, we evaluate semiclassically (q 2 oisc ) 
[23j], and insert it in the expression -E no isc = ^ K cs(q 2 i SC ) , 
with K e g the effective stiffness of the micro-mirror modi- 
fied by radiation pressure. Specifically, we first calculate 
the evolution of the operators describing the dynamics 
of the intracavity field and of the cantilever via a set 
of quantum Langevin equations (QLE) that include the 
hamiltonian dynamics of Eq. (|9|), input laser noise, cav- 
ity damping, and the Brownian noise associated with the 
coupling of the cantilever to a reservoir at temperature 
T 



a 

q = 
p 



1 

i(6 - gq)a - -a + ^ya m 7 
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Here 5 = lu c — uji is the detuning of the laser frequency 
toi from the cavity resonance w c , 7 is the cavity damp- 
ing rate and j c is the damping rate of the cantilever. 
The noise operator a m describes the laser field incident 
on the cavity. Its mean amplitude is (a m (t)) = a™ 
and its fluctuations are taken to be delta-correlated, 
(6a in (t)5a in ^(t')) = 5(t - t'). The Brownian noise opera- 
tor £ m describes the thermomechanical noise that couples 
to the cantilever from the environment. Its mean value 
is zero and, for high mechanical quality cantilevers such 
as considered here, ujq/j c ^> 1, its fluctuations are delta- 
correlated (8£ in (t)6C n (tj) = 7 c (2n + l)S(t - t'). 



The semiclassical, steady-state solutions of Eqs. fTO]) 
oscillate due to the time-dependent driving of the mir- 
ror. However the amplitude of these oscillations is small 
as compared to their mean value (for example with the 
parameters given in the caption of Fig. [21 the mean 
power inside the cavity is about 1.1 mW with an os- 
cillation amplitude of the order of lOOnW) about that 
value. We therefore assimilate the steady state solu- 
tions to their temporal mean values which are readily 
obtained as p s = and q s = g\a s \ 2 /ujq. The mean steady 
state field amplitude a s is complex in general, but it 
is possible to choose its phase such that a s is real, i.e. 



V7k n |/[(i) 2 + (<5-S<? s ) S 
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This amplitude 



can be bistable for high incident fields [2J] and appro- 
priate cavity-laser detunings, but one can easily choose 
system parameters that avoid this regime. 

To account for the quantum fluctuations 5a about the 
semiclassical steady-state a s we decompose a = a s + 5a, 
and similarly for the other operators of Eq. ([9]). Lin- 
earizing the quantum Langevin equations of motion in 
the fluctuations gives then 



v(t) = Av(t)+v(t), 



(11) 



where v(t) is the vector of fluctuations. Recasting the 
optical fields in terms of their quadratures X a and Yf, 
with fluctuations 5X a = (5a + 5a))/\2, 5Y a = (5a — 
6af)/iy/2], and similarly for the input field, the input 
noise vector becomes n(t) = (0, <5£ m , y/jSX" 1 , y/y5Y™) 
and 



A = 



LO 

-ojo -7c 



\ G 




G 

-7/2 A 
-A -7/2/ 
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where A = 5 — gq s is the effective cavity detuning and 
G — g s V% is the effective optomechanical parameter. We 
have verified that the Routh-Hurwitz criterion is satisfied 
here, that is, the eigenvalues of A all have a negative 
real part which ensures the stability of the steady-state 
solution. 

We can determine (q 2 oisc ) = h(q 2 ) / Mcu from the cor- 
relation matrix C, whose elements are given by = 
[(uj(oo)uj(oo) + Vj(ao)vi(ao))]/2. This matrix can be 
shown to be the solution of the equation AC + CA T = 
—D, with D = diag[0,7 c (2n+ 1), 7/2,7/2] [H]. The full 
solution of this equation can be obtained analytically, 
but it is cumbersome and will not be reproduced here. 
For our present purpose, we only need to consider the 
first diagonal element of the matrix C as it yields (q 2 )- 
From this, the integral of Eq.© has been evaluated an- 
alytically in the limit Q ^> 1. It is illustrated in Fig. [2] 
for the same parameters as in Fig. [T] except M = 50fg 
(l/xmx200nmxl50nm silicon cantilever), a cavity length 
of L ~lmm, ALascr = 810nm, a cavity resonance fre- 
quency lo c ~ 2tt x 10 14 Hz and effective cavity detuning 
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FIG. 2: (Color online) Cavity mass sensor mass sensitivities 
in Dalton units (lDa~ 1.7 x 10 _27 kg) as a function of temper- 
ature for the parameters of the text and P = 100/iW (dashed 
red line) and P = lmW (solid blue line) and T = 5. The dot- 
ted and dot-dashed lines are the mass sensitivities at T = OK 
respectively for P = 100/iW and P = lmW. 



A = ujq. The laser power P and the finesse J- are given in 
the caption. Note that with the further approximation 
27rA/ <C uo/Q and with Eq.© we obtain a simplified 
expression for the mass sensitivity 



SM S = 



E d 



1/2 



(An 
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For a cavity finesse T = 5, the level of noise intro- 
duced by the cavity setup equals the noise due to the 
vacuum fluctuations of the cantilever. The mass sensi- 
tivity is therefore reduced by a factor of two from that of 
the free cantilever mass sensor when it is a part of this 
optical arrangement. One attractive feature of the cavity 
mass sensor is the weak dependence of its sensitivity on 
temperature. Indeed, although increasing the input laser 
power increases the level of noise at T = OK as more light 
is present in the cavity, it decreases the temperature de- 
pendence of the mass sensitivity, so that for P = lmW 
the mass sensitivity at T = 25K is only about 600% 
worse than it is at T = OK. This is to be contrasted to 
the situation for the free cantilever mass sensor for which 
the mass sensitivity degrades by a factor of 30 over the 
same temperature range. This drastically different be- 
havior is a result of the active cooling by the laser which 
maintains the micro-mirror at an effective temperature 
much smaller than that of the bath 25 1. 

On the other hand, one drawback of the cavity mass 
sensor is the small mass required of the micro-mirror. 
Eq. ((3J) shows that the smallest detectable mass change is 
proportional to the mass of the cantilever. For our pa- 
rameters it follows that a cavity sensor mass of M ~ 0.1 fg 
is needed for SM S to approach an electron mass. We also 
find that the cavity mass sensor requires a cavity finesse 
approaching T ~ 5 to reach a regime in which the mass 
sensitivity is limited by vacuum fluctuations of the can- 
tilever, implying a reflectivity R > 30%. It might be dif- 
ficult to combine both characteristics - small mass and 
high reflectivity - in one cantilever and we have there- 



fore considered a 50fg, l^tmx200nmx 150nm silicon can- 
tilever for the cavity based mass sensor. Future work we 
will study whether a three-mirror configuration presents 
additional advantages [5^, [2(| . 

As a final remark, we note that more advanced de- 
tectors based on squeezed cantilever motion will only be 
attractive if the level of noise introduced by the optical 
setup that produces such nonclassical cantilever states 
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22| drops below the noise level due to vacuum fluc- 
tuations of a cantilever in a coherent state (that is for 
T > 5). An alternative might be to turn off the laser 
when the cantilever motion is squeezed. This would 
require stroboscopic mass measurements which would 
in turn compete with the requirements on the value 
of Af needed to obtain good mass sensitivities (here 
A/ ~lkHz). 

In summary, we have determined the quantum limit to 
nanocantilever-based mass sensing to be of the order of 
an electron mass for realizable parameters; we also note 
that the mass sensitivity of the quantum cantilever, as 
opposed to its classical counterpart, can be independent 
of its frequency at low temperatures. Our findings are 
of use to high sensitivity detection of isotopes and sto- 
ichiomctry of compounds. It will also be intriguing to 
explore the consequences of our findings for spin, charge 
and force sensing. 

This work is supported in part by the US Office of 
Naval Research, by the National Science Foundation and 
by the US Army Research Office. 
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